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ABSTRACT 

The  effect  of  seeding  shock  waves  with  alkali  metal  smoke  is  studied.  The 
pressure,  density,  temperature  ratio,  electron  density,  and  conductivity  were 
calculated  with  a  digital  computer  and  tabulated  for  various  combinations  of 
gas  and  seeding  material.  A  hydrogen-driven  shock  tube  was  constructed  to 
measure  the  properties  of  the  seeded  gas  behind  the  shock  wave.  The  conduc¬ 
tivity  of  argon  seeded  with  cesium  was  measured  with  a  conductivity  probe  and 
compared  with  theoretical  calculations. 

1.  INTRODUCTION 

The  interaction  of  ionized  gas  flow  with  a  magnetic  field  has  been  known  for 
some  time  and  considerable  effort  has  been  made  to  investigate  and  utilize 
this  phenomenon.  The  effect  of  magneto  gas  dynamics  has  been  applied  to 
power  generation,  nuclear  fusion,  and  space  propulsion.  The  magnitude  of  all 
such  interactions  depends  on  the  electrical  conductivity  of  the  gas. 

One  method  of  studying  this  interaction  is  to  use  the  high  temperature  and 
hence  ionized  gas  behind  a  shock  wave  in  a  shock  tube.  However,  all  gases  of 
engineering  importance  have  comparatively  high  ionization  potentials  and 
hence  only  small  conductivities  can  be  obtained  easily.  One  way  to  increase 
the  interaction  is  to  increase  the  conductivity  by  an  appropriate  rise  in  the 
shock  Mach  number.  This,  however,  reduces  the  interaction  time  and  creates 
experimental  problems.  An  alternative  way  is  to  increase  the  conductivity  by 
introducing  small  quantities  of  easily  ionizable  material  into  the  gas.  The 
alkali  metals  or  alkali  earth  groups  have  low  ionization  potentials  and  are, 
therefore,  suitable  as  seeding  material.  The  seeding  materials  are  easily  ionized 
in  the  hot  gas  behind  the  shock  wave  and  contribute  a  large  number  of  electrons, 
thus  increasing  the  conductivity. 

In  the  present  report,  the  effect  of  seeding  shock  waves  in  the  inert  gases 
argon  and  helium  with  alkali  metals  potassium  and  cesium  was  investigated. 
The  inert  gases  were  chosen  to  avoid  the  complexity  of  the  dissociation  problem. 

2.  SEEDED  SHOCK  WAVES 

The  gas  properties  behind  a  moving  normal  shock  wave  in  various  gases 
have  been  calculated  in  a  number  of  previous  investigations  (Glass  1958). 
For  a  monatomic  gas  which  is  seeded  with  small  quantities  of  alkali  metals,  the 
analysis  is  quite  similar.  As  a  simplification,  the  ratio  of  specific  heats,  y,  is 
taken  to  be  5/3  for  both  the  gas  and  the  alkaline  metal  vapor.  This  is  justified 
since  the  excitation  energy  is  only  a  small  portion  of  the  total  enthalpy  in  the 
present  range  of  work. 
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We  define  the  degree  of  seeding  by 

,n  w,  +  n,+  _  p, 

ng  +  «gT  pg  Ma  ’ 

where  n  is  the  number  density,  M  is  the  molecular  weight,  g  denotes  the  gas  and 
s  the  seed  material.  Assuming  an  equilibrium  state  (2)  behind  and  a  quiescent 
state  (1)  in  front  of  the  shock,  the  properties  of  regions  (1)  and  (2)  are  related 
by  the  conservation  of  mass,  momentum,  energy,  the  Saha  equation,  and  the 
equation  of  state: 

(2)  ft,  -  £  - 1  +  7 M (i  +  f;)(i  -  £) . 

(3)  Hi  +  |(w  -  «2)2  =  Hi  +  *«/2. 
where 
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=  !  (1  +  «g)  17-  T2  +  ctgjf-  +  «.-TT-c  +  5  (1  +  «■>)  7/  7 
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Mg 
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The  above  equations  have  been  solved  simultaneously  for  Pn,  921 ,  T2i,  and 
ne  on  a  computer  for  5  =  c.MJMg  =  0,0.0002,0.002,  and  0.02.  The  pressure, 
density,  and  temperature  differ  only  very  slightly  between  the  seeded  and  the 
unseeded  case.  The  electron  density  ne,  however,  is  increased  appreciably  and 
is  plotted  in  Fig.  1. 


3.  ELECTRICAL  CONDUCTIVITY  BEHIND  SHOCK  WAVES 

The  electrical  conductivity  of  shock-heated  argon  has  been  measured 
experimentally  and  computed  theoretically  by  Lin  et  al.  (1955).  There  the 
reciprocal  of  the  conductivities  due  to  close  encounters  (between  electron  and 
neutrals)  and  distant  encounters  (between  electron  and  ions)  was  added  to 
give  the  reciprocal  of  the  total  conductivity.  The  average  collision  cross 
section  was  used  instead  of  correctly  using  the  total  collision  frequency 

V total  (»)  =  "elW  +  Vm(v), 

and  electron-electron  collisions  are  neglected.  To  obtain  a  better  approxi¬ 
mation,  the  analysis  of  Shkarofsky  (1961)  was  used.  In  this  reference,  the 
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Fig.  1.  Electron  density  of  shock  waves  in  seeded  argon  and  helium. 


collision  frequencies  between  electrons  and  neutrals  were  assumed  to  vary  as  a 
simple  power  r  of  the  electron  velocity.  The  tensor  conductivity  including  all 
interactions  is  given  by 


(10) 


/  b  +  c  j(b  —  c) 

I  —j(b  -  c)  b  +  c 

\  0  0 


where 
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b  = 


rife1 _ 1 _ 

2 m,  ((I'm)  +  (>'el»g  +i(w  —  03b)h  ’ 
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«ee 
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m  ((vm)  +  (>'ei))g  +  joih 
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/oa  =  Maxwellian  distribution  function  of  electrons  and 
g  and  It  are  functions  of  ub/((vm)  +  <t-ei>),  (vm)/{ve{),  and  Z. 

,  v  4(2ir)1/2Meln  A  ( 

“  W'\kl r/2  \4 ireo/  ’ 

3  (47T60ifer)3/2 
a_2 z?  (x«;)t/2  • 


Z  =  charge  number  =  1  in  the  present  case. 


The  values  of  g  and  h  were  computed  and  plotted  as  functions  of  («  ±  ub)/ 
((I'm)  +  (fci))  and  (vm)/(i' ci)  in  Shkarofsky  (1961)  for  values  of  r  =  —3,  —2, 
—  1,  0,  1,  2,  3.  A  review  of  the  literature  on  momentum  transfer  cross  sections 
gives  the  approximate  dependence  of  the  collision  frequencies  between  neutrals 
and  electrons  on  the  electron  energy  (Frost  1961).  The  collision  frequencies 
for  the  various  seed  and  gas  materials  used  in  this  report  are  given  below. 
For  the  case  of  argon,  the  expression  is  derived  from  a  curve  of  the  momentum 
transfer  cross  section  obtained  by  Frost  and  Phelps. 

He  vm/N  =  3.14  «1/2 

K  vm/N  =  160 

Cs  vm/N  =  160 

A  Q  =  9.612X10-17  «'  35(0.3«Xl0eV) 

N  =  number  of  neutral  atoms. 


In  the  seeded  gas  behind  the  shock,  therefore,  the  collision  frequency  between 
electron  and  neutrals  is 

(11)  Fm  =  A  +  BVC. 

Since  the  results  presented  in  Shkarofsky  are  based  on  the  power  law  relation 
vm  =  constant  X  vT  (r  =  —3,  —2,  —  1, 0,  1,  2,  3),  it  is  necessary  to  approximate 
the  collision  frequency  by  a  similar  expression  in  order  to  use  the  tabulated 
results  in  g  and  h.  Assuming  that  the  relation  is 

(12)  Vm  =  AV, 

the  values  of  the  constants  K  and  r  can  be  obtained  in  terms  of  .1 .  B,  C,  and  t'o 
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if  we  take  only  small  variations  of  v  from  the  center  point  v0.  v0  is  then  the 
electron  velocity  corresponding  to  the  temperature  at  which  we  want  to  com¬ 
pute  the  conductivity.  Equating  equations  (11)  and  (12),  and  assuming 
v  =  t>0(l  +  £)  where  £  «  1,  vn  is  approximated  by  the  equation: 

=  (A  +  BvoC )  (  -  )  BCv0C/(A  +  Bv0C), 

\  Vo  / 


and 


(vm) 


4_7T 
3  ne 


vmdv 


feT(2.5  +  \r)  ( 2k]\n 

r(2.5)  V  m  )  ' 


Using  the  above  equations,  (pe))  and  (vm)  were  calculated  for  each  Mach  number 
and  gas-seed  combination  on  a  digital  computer.  The  value  of  v0  depends  on  the 
equilibrium  temperature  and  is  calculated  as  follows.  The  d-c.  conductivity 
without  magnetic  field  and  considering  neutral  scattering  is  given  by 


4re2  pi  3df0, 

<r  - z —  I  —v  ~r~  dm 

3m  Jo  vm  dv 

To  find  the  power  law  for  vm  that  would  best  suit  the  gas  and  seed  mixture, 
the  center  point  t»0  is  taken  to  be  at  the  maximum  of  the  integrand  (v3(df0/dv)). 
Differentiating  the  integrand  and  setting  it  equal  to  zero,  and  assuming  that 
/0  is  Maxwellian, 


4vV 


mv  /2kT 


m  5  —mA /2k  T 

kfv  e 


=  o, 


from  which  i>02  =  4 kT/m.  The  values  of  g  were  interpolated  from  the  results  in 
(2)  and  a  is  calculated  by  the  following  equation: 

2  , 

_  WeC _ 1 _ 

m  g({vm)  +  qvi»  ' 

The  effect  of  a  large  magnetic  field  can  be  incorporated  without  difficulty.  The 
values  of  g  and  h  can  be  interpolated  as  above  and  the  tensor  conductivity 
calculated  according  to  equation  (10). 


4.  EXPERIMENTAL  APPARATUS 

The  shock  tube  used  in  this  investigation  consists  of  a  cylindrical  driven 
section  about  30  feet  long  and  a  driver  section  about  G  feet  long.  Both  sections 
have  inside  diameters  of  2|  inches.  The  test  section  is  a  cylindrical  pyrex  pipe 
about  2  feet  long  with  inside  diameter  of  I5  inches.  Ahead  of  the  pyrex  pipe 
is  a  4-foot,  1^-inch  inside  diameter  steel  section  which  can  be  inserted  into  the 
driven  section  of  the  shock  tube.  This  has  the  effect  of  discarding  the  !  urbulent 
boundary  layer  induced  behind  the  shock,  together  with  any  di  iphragm 
fragments  or  impurities  that  might  have  been  swept  from  the  shock  tube  wall. 
The  insert  prevents  the  reflected  shock  wave  from  disturbing  the  flow  into  tin 
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test  section.  A  dump  tank  is  connected  to  the  end  of  the  test  section  to  relieve 
the  pressure  and  prevent  the  formation  of  a  reflected  shock.  A  metal  plunger 
activated  by  compressed  air  was  designed  to  break  the  diaphragm.  Several 
types  of  metal  diaphragms  were  tested  and  copper  diaphragms  up  to  0.0:52  inch 
thick  were  found  to  be  best  since  they  petal  nicely  and  no  fragments  were 
broken  off.  The  experimental  setup  is  shown  in  Fig.  2. 

With  driver  pressure  up  to  1000  p.s.i.  hydrogen  and  channel  pressure  from 
2.5  to  20  mm  Hg,  shock  Mach  numbers  from  6  to  10.5  were  obtained  in  the  test 
section.  Nitrogen  was  used  to  flush  out  the  hydrogen  after  each  run. 

The  conductivity  of  the  hot  gas  behind  the  shock  was  measured  with  the 
conductivity  probe,  the  theory  of  which  is  described  in  detail  by  Savic  (19(52), 
and  only  a  brief  description  is  given  here.  The  probe  consists  of  a  (1-inch-long 
pyrex  tube  which  fits  over  the  test  section,  and  has  six  coats  of  platinum  paint 
baked  on  the  inside  to  eliminate  the  capacitive  coupling  effect  of  any  substance 
put  inside  the  tube.  Two  turns  of  copper  wire  were  wound  on  the  middle  of  the 
tube,  constituting  part  of  the  inductance  of  a  tuned  circuit.  This  inductance 
changes  slightly  when  the  hot  conducting  gas  behind  the  shock  passes  through 
the  test  section,  and  appears  as  a  frequency  change  in  the  tuned  circuit.  It  is 
then  fed  into  a  discriminator  circuit  which  transforms  it  into  a  d-c.  voltage. 
Mercury  and  sulphuric  acid  were  used  for  its  calibration  since  there  are  no 
substances  within  the  right  range  of  conductivity  except  for  fused  salts  and 
semiconductors.  For  this  reason,  considerable  difficulty  arises  in  obtaining  a 
reliable  calibration  and  the  accuracy  of  the  probe  is  limited  to  about  15%. 

5.  SEEDING  OF  THE  SHOCK  TUBE 

There  are  two  different  methods  of  seeding  the  shock  tube.  The  first  one  is  to 
inject  the  vapor  of  the  alkali  metals  into  the  gas.  This  is  the  ideal  way  of  seeding 
since  a  genuine  gas-seed  mixture  can  be  obtained.  However,  this  necessitates 
heating  the  shock  tube  (at  least  the  driven  section)  up  to  a  temperature  where 
the  vapor  pressure  of  the  alkali  metal  is  sufficiently  high.  An  oven  would  have 
to  be  constructed  to  enclose  the  whole  shock  tube  and  a  predetermined  amount 
of  alkali  metal  introduced  inside  the  tube.  This  raises  problems  involving 
measurements  of  the  shock  properties,  especially  optical  measurements  and 
material  problems,  since  alkali  metal  vapor  at  this  temperature  is  very  corrosive. 

An  alternative  method  of  seeding  the  shock  tube  is  by  introducing  a  cloud 
of  small  particles  of  alkali  compound  into  the  shock  tube.  This  can  be  done  In- 
blowing  the  gas  through  a  finely  divided  powder  of  the  seeding  material 
(Basu  1960)  or  bubbling  it  through  the  molten  material  and  allowing  the 
vapor  to  cool.  However,  the  above  methods  do  not  produce  a  cloud  ol  fine 
uniform-size  particles.  A  slight  improvement  is  obtained  by  heating  the  cesium 
chloride  to  about  1000°  C  and  letting  a  jet  of  cool  argon  blow  onto  the  vapor. 
The  compound  sublimes  readily  at  this  temperature  and  above  the  molten 
cesium  chloride  there  is  formed  a  region  of  high  concentration  of  small  crystal¬ 
line  particles.  The  argon  jet  blowing  on  the  surface  then  carries  oft  these 
particles,  which  are  suspended  in  the  form  of  smoke.  A  sample  of  this  smoke  was 
deposited  on  a  small  carbon-coated  grid  and  examined  under  an  electron 
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microscope.  It  was  seen  to  consist  of  minute  particles  having  linear  dimensions 
of  about  500  A,  together  with  large  particles  which  are  aggregates  of  a  few 
tens  or  even  hundreds  of  the  basic  crystals.  A  few  plates  with  various  magni¬ 
fication  are  shown  in  Fig.  5. 

The  shock  tube  and  the  seeding  apparatus  are  evacuated  initially.  The  CsCl 
in  a  quartz  container  is  then  heated  up  to  about  1000°  C.  A  jet  of  argon  is 
directed  onto  the  surface  of  the  molten  compound  and  a  predetermined  amount 
of  the  smoke  is  let  off  from  the  top  of  the  quartz  container  into  a  cylindrical 
chamber  which  is  separated  from  the  shock  tube  by  a  butterfly  valve.  A  sample 
bottle  with  two  electrodes  is  connected  to  the  cylindrical  chamber  for  measuring 
the  amount  of  seeding.  After  the  required  amount  of  smoke  has  been  fed  into 
the  chamber,  the  butterfly  valve  is  flipped  open  and  the  smoke  fills  the  shock 
tube  up  to  the  desired  pressure.  The  shock  tube  is  then  fired  immediately  to 
avoid  settling  of  the  smoke  particles  on  the  shock-tube  walls. 

The  degree  of  seeding  is  measured  by  adding  water  to  the  sampling  bottle 
and  measuring  the  conductivity  of  the  solute  between  the  electrodes. 

The  disadvantages  of  seeding  with  solid  smoke  particles  of  a  compound  are: 

1.  The  seeding  material  has  to  be  vaporized  and  dissociated  before  yielding 
electrons.  However,  in  the  present  range  of  work,  the  energy  required  to  vaporize 
and  dissociate  the  alkaline  compound  is  but  a  few  percent  of  the  total  enthalpy 
of  the  hot  gas,  so  that  the  effect  on  the  final  conductivity  is  very  slight. 

2.  The  affinity  of  the  nonmetallic  radical  on  electrons  leads  to  a  reduction  of 
the  electron  density.  However,  an  estimate  shows  that  the  reduction  of  electron 
density  bv  the  nonmetallic  radical  is  only  appreciable  at  temperatures  below 
5000°  K. 


6.  RESULTS 

The  theoretical  conductivities  of  seeded  inert  gases  are  plotted  in  Figs.  3  and 
4.  The  first  is  a  comparison  between  Lin’s  method  and  that  of  Shkarofsky. 
The  curves  tend  to  converge  at  both  ends  and  the  discrepancy  is  worse  at  the 
intermediate  range.  This  is  expected  since  Lin’s  analysis  sums  the  resistivity 
from  electron-neutral  and  electron-ion  encounters  which  are  accurate  at  the 
extreme  ends  of  slightly  ionized  and  fully  ionized  gases  respectively.  At  both 
ends,  one  of  the  two  is  negligible,  so  that  the  discrepancy  is  slight;  but  at  the 
intermediate  range,  the  two  types  of  encounter  cannot  be  summed  as  in  Lin’s 
analysis  since  their  velocity  dependence  is  vastly  different.  For  example,  in  the 
case  of  argon,  the  collision  frequency  is  proportional  to  r3  for  the  electron- 
neutral  encounters  and  to  v~3  (or  v~*  in  Frost,  taking  into  account  the  electron- 
electron  interactions)  for  the  electron-ion  encounters. 

From  the  theoretical  curves,  it  can  be  seen  that  the  seeded  shocks  have 
approximately  constant  conductivity  at  lower  shock  speeds,  and  the  con¬ 
ductivity  rises  to  meet  the  unseeded  case  at  higher  speeds.  This  is  because  at 
lower  .1/ft,  most  of  the  conductivity  is  contributed  by  the  seeding  material, 
which  is  practically  all  ionized  above  d/»  =  6,  whereas  the  gas  is  only  slightly 
ionized.  As  d/»  increases,  however,  the  gas  becomes  more  and  more  ionized 
until  the  percentage  ionization  exceeds  the  percentage  seeding,  and  from  then 
on  the  contribution  from  the  gas  predominates. 
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The  experimental  results  are  also  plotted  in  Figs.  4(a)  and  4 (b).  The  amount 
of  seeding  was  approximately  2%  since  the  actual  amount  inside  the  shook  tube 
cannot  be  measured.  However,  variation  of  the  amount  of  seeding  slighth 
above  or  below  2%,  does  not  produce  appreciable  difference  in  the  measured 
conductivity. 

The  measured  conductivity  agrees  fairly  well  with  the  theoretic  al  value  for 
the  unseeded  case,  whereas  for  the  seeded  shocks,  the  measured  values  arc-  to  > 
small  by  approximately  half.  This  may  be  due  to  a  number  of  factors.  Streak 
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Fig.  4(a).  Conductivity  of  argon  seeded  with  cesium  (pi  =  10  mm  Hg).  Fig.  4(b).  Conductivity  of  argon  seeded  with  cesium  (pi  —  5  mm  Hg). 


I  k;.  Me).  Conductivity  of  argon  seeded  with  potassium  (pi  =  10  mm  Hg).  Fig.  4(d).  Conductivity  of  argon  seeded  with  potassium  (pi  =  5  mm  Hg). 
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Fig.  4(e).  Conductivity  of  helium  seeded  with  cesium  (pi  =  10  mm  Hg).  Fig.  4(f).  Conductivity  of  helium  seeded  with  cesium  (pi  —  5  mm  Hg). 
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>'  Of  helium  seeded  with  potassium  {pi  =  10  mm  Hg).  Fig.  4(A).  Conductivity  of  helium  seeded  with  potassium  {pi  =  5  mm  Hg). 
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photographs  by  the  method  of  Petschek  and  Byron  (1957)  show  that  while  the 
relaxation  time  ranges  from  about  100  to  400  microseconds  for  the  unseeded 
shocks,  it  was  very  small  for  the  seeded  ones  (Fig.  6).  This,  together  with  the 
steep  rise  in  the  conductivity  trace  (Fig.  7),  suggests  that  the  low  measured 
conductivity  is  probably  not  due  to  nonequilibrium. 
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One  probable  cause  of  the  low  observed  conductivity  is  the  particle-size 
distribution.  The  largest  aggregates  have  a  settling  rate  of  about  La  cm/sec 
compared  with  about  3  X  10~3  cm/sec  for  the  small  particles,  so  that  at  the 
bottom  of  the  shock  tube  there  is  a  high  concentration  of  large  aggregates.  Even 
though  the  estimated  time  required  to  vaporize  and  dissociate  the  small 
particles  is  of  the  order  of  microseconds,  the  time  elapsed  before  the  large 
aggregates  can  contribute  to  the  conductivity  might  be  much  larger,  especially 
when  the  largest  ones  lie  in  the  cool  boundary  layer.  A  few  shots  were  taken 
with  deliberate  time  delay  between  seeding  and  shooting,  and  the  observed 
seeded  conductivity  is  even  lower  than  that  of  the  unseeded  case.  This  shows 
that  the  additional  energy  required  to  vaporize  and  dissociate  the  rich  aggregate 
of  particles  at  the  bottom  of  the  tube  tends  to  cool  down  the  boundary  layer 
even  further. 

The  diffusion  time  required  for  the  hoi  seeded  gas  to  form  a  homogeneous 
conducting  fluid,  assuming  ambipolar  diffusion,  ranges  from  a  fraction  to  about 
50  microseconds  for  particle  separation  of  100  and  1000  microns  respectively, 
so  that  this  effect  is  unlikely  to  cause  lower  observed  conductivities.  However, 
since  the  smoke  particles  can  only  be  examined  after  being  deposited  on  fine 
grid,  there  is  no  way  of  knowing  for  certain  what  the  separation  is  when  the 
particles  are  still  suspended  in  the  argon,  especially  when  they  tend  to  aggregate 
on  settling. 


7.  CONCLUSION 

The  theoretical  conductivity  of  argon  seeded  with  solid  smoke  particles  was 
not  fully  attained.  Because  of  the  uncertainties  of  using  solid  particles  as  dis¬ 
cussed  in  Section  6,  it  is  not  known  which  is  the  most  important  factor  in  the 
discrepancy  between  the  theoretical  and  measured  conductivities.  The  boundary 
layer,  the  solid  particle  distribution,  and  their  attainment  of  equilibrium 
ionization  can  be  investigated  with  optical  methods  such  as  schlieren  and 
spectroscopy  to  give  a  clearer  understanding  of  the  ionization  process  behind 
the  shock  wave.  However,  conductivities  of  argon  seeded  with  alkaline  salt  will 
be  measured  shortly  in  a  plasma  jet,  where  there  are  no  boundary-layer  and 
solid-particle  problems  since  the  salt  is  injected  as  a  vapor.  These  results  will 
then  confirm  whether  or  not  the  low  conductivity  is  due  to  the  effects  mentioned 
above. 

In  view  of  the  present  experimental  results,  it  can  be  concluded  that  even 
though  particulate  seeding  in  shock  tubes  does  not  give  the  desired  conductivity, 
it  does  raise  the  conductivity  appreciably  at  low  Mach  numbers.  The  main 
advantage  is  the  simplicity  of  adding  the  seed  material,  and  all  the  incon¬ 
veniences  of  heating  the  shock  tube  in  an  oven  are  avoided. 
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